REPORT 1158 


PREDICTION OF FLAME VELOCITIES OF HYDROCARBON FLAMES 1 

Bv Gordon L. Dugger and Dorothy M. Simon 


SUMMARY 

The laminar-flame-velocily data previously reported by the 
Lewis laboratory are surveyed with respect to the correspondence 
between experimental flame velocities and values predicted by 
semilheorelical and empirical methods. The combustible mix- 
ture variables covered are hydrocarbon structure (56 hydro- 
carbons), equivalence ratio of fuel-air mixture, mole fraction 
of oxygen in the primary oxygen-nitrogen mixture ( 0.17 to 0.50), 
and initial mixture temperature (200° to 615° K). The semi- 
tlieoretical methods of prediction considered are based on three 
approximate theoretical equations for flame velocity: the 
Semenov equation (bimolecular) , which is primarily based on 
conductive heat transfer between the flame and the reactants; 
the Tanford-Pease equation, which is based on the diffusion of 
chain carriers of the oxidation reaction into the reactants; and 
the Manson equation, which is a modification of the momentum- 
pressure-drop equation that does not include chemical kinetics. 
In each equation a semiempirical factor is used to bring the 
predicted values for a given variable and fuel into the best aver- 
age agreement with the data, so that the variation in the relative 
prediction of the individual datum points may be considered. 
For the resulting semitheorelical equations, it is assumed that 
thermal equilibnum is attained at the end of the flame zone, and 
values for the transport properties are estimated by extrapolation 
and simple additive relations. The empirical relations between 
combustible mixture variables and flame velocity are based on 
the usual methods of correlation. 

Within these limitations, the results may be generalized 
as follows: 

(1) The three semitheoretical equations predict relative 
flame velocities reasonably well, generally with mean 
deviations of 2 to 15 percent. 

(2) Considering only the bimolecular fuel-oxygen reaction 
with low-temperature activation energies, the Semenov 
equation may be used to give good relative predictions when 
a “steric factor” is determined semiempirically. 

(3) The average “specific rate constants” k a obtained 
with the Tanford-Pease equation for molecular fuel — active 
particle reactions behave as Arrhenius rate constants only 
with the data orf the effect of initial temperature, not with 
the equivalence-ratio or oxygen-concentration data. The 


relative predictions obtained by this equation are not very 
sensitive either to the temperature dependence assigned to 
the diffusion coefficients or to the recombination factor 
computed for hydrogen atoms. Regarding active particles 
considered, better results are generally obtained for hydro- 
carbon flames when the hydroxyl radical and the hydrogen 
and oxygen atoms (OH, H, and O) are considered rather 
then H alone; this is particularly true for studies of the 
effect of equivalence ratio. 

(4) For the equivalence-ratio and oxygen-concentration 
data, considerably better relative predictions are obtained 
from the Manson equation when the pressure drop across 
the flame front is considered to be due to H, OH, and O 
than when only H atoms are taken into account. 

(5) The example presented by the ethylene data shows 
that an empirically determined rate constant or proportion- 
ality factor from equivalence-ratio data at atmospheric 
pressure and room temperature may predict the effect of 
initial temperature or oxygen concentration within approxi- 
mately 20 percent by the Semenov or Tanford-Pease equation 
and within approximately 30 percent by the Manson 
equation. 

(6) For engineering applications, the effects of the 
parameters studied could be estimated’just as satisfactorily, 
and more easily, by one or another of the empirical correla- 
tions indicated, as compared with the three semitheoretical 
equations considered. However, the use of the semitheo- 
retical equations in some cases reduces the number of 
constants required. 

INTRODUCTION 

The ability to predict flame velocities of fuels is of growing 
importance in the field of aircraft propulsion, since a cor- 
relation has been found between the combustion efficiency 
of a ram-jet burner and the laminar flame velocity of the 
fuel (ref. 1). The prediction of flame velocities is difficult 
for three reasons: (1) There is no complete, rigorous theory 
which can be readily applied. There are, however, a number 
of approximate equations in the literature which approach 
the problem of flame propagation from various viewpoints. 
(2) There are no data on the kinetics of the oxidation process 
under flame conditions and very few data on transport 


1 Supersedes NAOA RM E52J13, “Prediction of Flame Velocities of Hydrocarbon Flames," by Gordon ~L. Dugger and Dorothy M. Simon, 1952. 
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properties at liigli temperatures. (3) Different methods of 
flame-veloci ty measurement give different values; thus, it 
is difficult to compare data from different sources. The 
uncertainty in measurements made by a given method is of 
the order of 5 percent. 

The present stud} 7 and the calculations required for it 
were made during 1951-52. The study .covers all the data 
obtained on the effects of hydrocarbon structure, initial 
temperature, and hydrocarbon-oxygen-nitrogen mixtures 
which were obtained at the NACA Lems laboratory prior 
to September, 1952 (refs. 2 to 9). The semitheoretical 
methods are based on the Semenov equation (ref. 10), the 
Tanford-Pease square-root law (refs. 11 and 12), and the 
Manson equation (ref. 13). These three equations were 
derived from different models of the flame-propagation 
process. The Semenov model is essentially a thermal model 
which includes chemical reaction kinetics; the Tanford- 
Pease model is based on the diffusion of chain carriers of 
the oxidation reaction; and Manson used a modification of 
the momentum-pressure-drop equation which does not in- 
clude chemical kinetics. The empirical relations for the 
effect on maximum flame speed of hydrocarbon structure, of 
initial mixture temperature, and of oxygen concentration 
were based on the usual methods of correlation. 

It is recognized that there are other approximate theoreti- 
cal equations which may give as good or better predictions 
of the NACA experimental results. The particular equations 
used were chosen because they exemplify three different 
approaches to the problem and because the Semenov and 
Tanford-Pease equations were used in previous NACA 
papei's and many of the calculations had already been made 
individually. It is also recognized that there are many 
flame-speed data in the literature which could be used in 
such studies; the present paper is confined to NACA data 
because many of the calculations had been made, and because 
it was desirable to avoid the complications arising from 
comparing data obtained by man} 7 different methods and 
teclmiques. 

NOMENCLATURE 

The following nomenclature is used in this paper: 
a fuel concentration, molecules of fuel per cm 3 of 

mixture 

b oxygen concentration, molecules of oxygen per 

cm 3 of mixture 

Bi term near unity arising from radical recombination 

C p molar heat capacity at constant pressure, cal/ 

(mole) (°K) 

c p specific heat, cal/ (g) (°K) 

c p mean specific heat, T 0 to T f , cal/(g) (°K) 

I) diffusion coefficient, cm 2 /sec 

Dt diffusion coefficient of itb active species into un- 

burned gas, cm 2 /sec 

E activation energy, kcal/g-mole 

E q fitted activation energy for group of data points, 

kcal/g-mole 

K empirical proportionality constant between ex- 

perimental flame velocity and value predicted 
by Manson equation 


K n empirical K from Manson equation when 
Ap=lJ2p E 7^r 

K-s empirical K from Manson equation when 

AP=l/2 (ph+Poh-^+Po-^J) 7jT 

K m Kz average K n or ify for group of data points 
k specific rate constant, cm 3 / (molecule) (sec) 

k a weighted mean k t for the three active particles II, 

OH, and O, each, reacting with fuel molecules, 
cm 3 / (molecule) (sec) 

ic a - average k a for group of data points, cm 3 / (molecule) 
(sec) 

k f specific rate constant for reaction between fuel 

molecule and ith active particle, cm 3 / (molecule) 
(sec) 

ki average k t for group of data points, cm 3 / (molecule) 

(sec) 

L m total concentration of gas at mean combustion- 

zone temperature, molecules/cm 3 
M molecular weight 

n total number of molecules of H 2 0 and C0 2 in prod- 

ucts per molecule of fuel by stoichiometric 
relation 

njno moles of reactants per mole of products from 
stoichiometric equation 

P steric factor, or probability factor, from the ex- 

pression k—PZ exp [ —E/RT] 

P average P for group of data points 

p total pressure of mixture 

Pi - mole fraction of ith active particle in burned gas 
(equilibrium flame temperature) 

Q mole fraction of potential combustion product in 

unburned gas 

Q! mole fraction of fuel in unburned gas 

R universal gas constant, kcal/(g-mole) (°K) or 

ergs/(g-mole) (°K) 

r exponent showing molecularity of flame reaction 

T absolute temperature, °K 

U flame velocity, cm/sec; implies maximum flame 

velocity -with respect to equivalence ratio, except 
in discussion of variation of equivalence ratio 
w reaction rate, molecules reacting per cm 3 /sec 

Xj mole fraction of jth component 

Z collision number; number of molecular collisions 

per second when concentration is one molecule 
of each type per cm 3 

a mole fraction of oxygen in oxygen-nitrogen portion 

of mixture, 0 2 /(0 2 -(-N 2 ) 

S percentage mean deviation in ratio of predicted 

flame velocity to experimental value for given 
group of data 

i) viscosity of mixture, poise 

rjj viscosity of jth component, poise 

6 m ratio of mean reaction-zone temperature to initial 

temperature 

X thermal conductivity, cal/ (cm 2 ) (sec) (°K/cm) 

p density of mixture, g/cm 3 
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a collision diameter, cm 

<p equivalence ratio, fraction of stoichiometric fuel- 

oxygen ratio 

Subscripts: 

a average value for the three active particles (H, 

OH, and 0) 

ejf effective value at mean reaction temperature 

/ conditions at computed equilibrium flame tempera- 

ture 

j any component of reaction products 

m conditions at mean reaction-zone temperature 

0 initial conditions 


APPROXIMATE THEORETICAL EQUATIONS USED 

SEMENOV BIMOLECULAIt EQUATION 

Zeldovich and Frank-Kamenetsky obtained an approxi- 
mate solution for the rate of flame propagation from the dif- 
ferential equations for heat conduction and fuel concentra- 
tion change across the flame front. This approximate solu- 
tion was presented in detail by Semenov (ref. 10). The 
general approximate solution may be written (see ref. 10, 
pp. 31, 48, and 49): 


V a 0 paC P (T f 


wdT 

-To) ( 2 ) 


( 1 ) 


where, for a bimolecular reaction between unlike molecules 
r=2 and 


j* f w dT=j* f k(a)(b) dT 


J>rn 2 

=a e f/b ef/ PZ -jjf- exp (— E/RT, ) 


( 2 ) 


For lean mixtures, ?><1, a e ti and b cff are computed by 
(following the assumptions of ref. 10, p. 45): 


a ef/= 

S* ' 

O 

II 

.(3a) 

h A To 

°eff — C 0 jT 

1 —<p 1 


(3b) 

For rich mixtures, <p~)> 1: 




T 0 

Ue// tto rj- 

1— - 
L 9 

o-ff)] 

(4a) 

be/r 

, To RT/IE . 
- a T f T f -T 0 

(4b) 


For stoichiometric mixtures, equations (3) and (4) are 
identical. All of the factors in equations (1) through (4) 
can, with the exceptions of P and E, be estimated by the 
extrapolation of thermodynamic tables (refs. 14 and 15) or 
calculated by the relations 

v—IDvjXj (5a) 


(V alues of rj calculated for combustion-product mixtures by 
equation (5a) were within 1 percent of values calculated by 
the method of ref. 16.) 


Op 

j 

(5b) 

X=(<7 p +|e)^= (ref. 15) 

(5c) 

22=1.336 — (ref. 17) 

p 

■ (5d) 

pM 
p ~~ RT 

(5e) 


(5f) 


Flame-velocity predictions by the Semenov equation are 
evaluated herein as follows. Activation energies from low- 
temperature reactions and measured flame velocities are 
substituted in equation (1) and a steric factor P is calculated 
for each experimental flame velocity. These P’s are then 
averaged to give P for the group of data points. The ratio 
of the predicted flame velocity to the experimental velocity 
is calculated as ( PjP ) 1/2 . The average deviation of the 
ratios from 1 is considered to be a measure of the accuracy of 
the predicted flame velocities. (The word accuracy is used 
to denote the degree of correspondence between flame ve- 
locities calculated by the methods described and measured 
flame velocities.) 

The data for various hydrocarbons consisted of maximum 
flame velocities (maximum with respect to equivalence ratio) 
of 56 pure hydrocarbons at 298° K and atmospheric pressure 
(refs. 2 and 3). These data were obtained by the NACA 
tube method. The hydrocarbons included straight and 
branched chain members of the alkane, alkene, alkadiene, 
alkyne, and cycloalkane series and benzene. No calcula- 
tions are made by the Semenov equation for these data, be- 
cause the work required by this treatment would not be 
warranted in view of the slight differences in flame velocities 
observed for most of the S ydrocarbons studied. However, 
it has been shown by more approximate calculations (ref. 18) 
that relative predictions would be satisfactory. Examples 
of these calculated values are given in table I. 

The initial-mixture-temperature data included maximum 
flame velocities of methane-air, propane-air, and ethylene-air 
mixtures obtained by a Bunsen burner method (total-area 
method, outside edge of shadow cast by flame cone, ref. 5) 
at seven initial mixture temperatures ranging from 200° to 
617° K. The following activation energies, in kilocalories 
per gram-mole, are reported in references 19, 20, and 21, 
respectively: methane, 51; propane, 38; ethylene, 40. These 
values are used in the Semenov equation to evaluate the 
flame-velocity predictions 
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TABLE I.— EXAMPLES OF PREDICTED MAXIMUM FLAME 
VELOCITIES FOR VARIOUS HYDROCARBON-AIR MIX- 
TURES AT ROOM CONDITIONS 



Flame velocity, cm/sec 

Hydrocarbon 

Measured 
(ref. 3) 

Calculated by 
simplified 
Semenov equa- 
tion (ref. 18) 

Calculated by 
Tanford-Pease 
equation using 
k a (ref. 25) 


33.8 

40.2 

33.9 


39.0 

39.0 

37.7 


38.5 

37.2 

37.3 


36.6 

38.8 

37.1 


36.3 

39.1 

37.0 


35.9 

37.0 

36.6 


43.8 

48.9 

45.3 


42.1 

42.0 

41.7 


39.3 

41.7 

40.7 


69.9 

64.9 

61.3 


48.5 

46.0 

48.6 


38.7 

38.3 

37.4 


40.7 

44.5 

40.4 







Figure 1. — Comparison of measured and calculated burning velocities 
for ethylene-air mixtures at different initial temperatures. 
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Type of data 

Derivation of empirical correlation, 
percent 

Deviations of semitheoretical predictions, percent 

Form of correlation 

Deviation 

Semenov 

bimolecular 

equation 

Tanford-Pease equation 

Manson equation 

Using F n 

Using k a 

Using K 
from Tf, 
data c 

Using K a 

Using K s 

Hydrocarbon structure ft 

U=NaKa+N b Kb+ . . . 

1.9 


3.3 

3.5 


3.6 

3.8 

Initial temperature b (200° to G15° K) 
Methane 
Propane 
Ethylene 

U=a+bT 0 n 

2.1 

1.6 

1.4 

2.7 
2.0 

1.8 

11.8 

12.5 

9.4 

7.8 
6.1 

5.8 

0.8 

1.8 

10.3 

10.0 

12.0 

10.1 

13.0 

13.8 

Equivalence ratio (range limited, see text) 
Methano, 307° K h 
Pentane * 

Propane, 302° K b 
Ethylene * 

U=a(Iog B)+6 . 

4.9 

2.2 

7.4 

3.0 

4.0 

5.0 

25.8 

64.8 

16.6 

5.7 

3.8 
0.6 
3.7 

7.0 

3.3 

23.6 

39.9 

15.1 

10.0 

1.8 

3.6 

1.2 

1.8 

Oxygen concentration » 

Propane, 311" K 
422° K 
Ethylene, 311° K 
422° K 

2,2,4-Trlmethylpcntanc, 311° K 
422° K 

|u=oT 0 "(a-6) 

J U=aT<?( a -b) 

ju=a(cc-6) 

U*=aTf(*-b) 

5.5 
3.2 

J 2.0 

l 2.7 

2.6 

{ 17.3 

1 12.4 

{ 4.0 

1 3.0 

6.4 
7.3 


6.4 
5.2 
4.9 
4.6 

4.5 
4.0 


16.5 

12.7 

7.5 

6.8 

7.3 

8.2 

15.2 

13.0 

5.2 

5.3 
1.1 
1.1 


» Data determined by tube method at room temperature (refs. 2 to 4). 
b Data determined by Bunsen burner method, shadow cone (ref. 5). 

0 Data determined by Bunsen burner method, schliercn cone (rets. 0 and 7). 


An example of these relative predictions is presented for 
ethylene-air mixtures in figure 1. Using the average steric 
factor P for this group of data and the aforementioned acti- 
vation energy results in perfect agreement between the middle 
portion of the curve predicted by the Semenov bimolecular 
equation and the experimental points. The predicted curve 
diverges slightly from the experimental points at higher or 
lower flame velocities so that the mean deviation for the 
whole range is 1.8 percent. The average deviations obtained 
for the three fuels are tabulated in table II under “Semenov 
equations.” It may be noted that essentially the same aver- 
age deviations would result from the curves in reference 5 
although the calculations were made by a more approximate 
method and the predicted flame velocities were referred to 
the experimental value for 25° C. 

Two kinds of mixtm-e-composition variables are studied. 
The first is the change of hydrocarbon concentration in air 
over an equivalence-ratio range of 0.7 to 1.3. Flame-velocity 
values by the tube method for ethylene-air and pentane-air 
mixtures at 298° K (ref. 4) and Bunsen-burner values (total- 
area method, outer edge of cone shadow) for methane-air mix- 
tures at 307° K and propane-air mixtures at 302° K (ref. 5) are 
used. When the effect of changes in equivalence ratio on the 
flame velocity is computed by the Semenov equation, it is 
found that the approximations used in equations (3) and (4) 
are not consistent for the region near stoichiometric, partic- 
ularly for equivalence 'ratios between 0.95 and 1.05. This 
fact is shown by the difference between the Semenov curves, 
in figure 2 (a), where points in the range 0.95 to 1.05 are 
included in plotting the curve, and figure 2 (b), where these 
points are omitted. If this critical region is omitted, the 
accuracy of the predicted flame velocities over the equivalence 
range is 3 to 7 percent. Specific deviations are listed in 
table II. 


The second composition variable studied is the ratio of 
oxygen to oxygen plus nitrogen a=0 2 /(0 2 +N 2 ). Maximum 
flame velocities (Bunsen-burner total-area method, outside 
edge of schlieren cone, refs. 6 and 7) for several a ratios were 
available for the following fuels and ranges of a: propane, 
0.17 to 0.50; ethylene, 0.17 to 0.35; and 2,2,4-trimethyl-” 
pentane, 0.21 to 0.50. The predicted curves from references 
6 and 7 are reproduced in figure 3 (a). As in the initial- 
temperature example in figure 1, the predicted curves cross 
the experimental curve near the middle of the flame-velocity 
range. The accuracy of prediction of these data is 3 to 17 
percent (table II). 



(a) Semenov bimolecular equation (thermal mechanism). 

Figure 3. — Effect of oxygen concentration on maximum flame velocities. 
Comparison of theoretically predicted curves with experimental results. 
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It may be noted in table II that flame-velocity predictions, 
by tiie Semenov equation for ethylene flames are accurate to 
5 percent for variations in both composition and initial 
temperature. For these predictions (as discussed pre- 
viously), each type of data was considered separate!} 7 , and 
sterie factors were calculated from the data. Now usually, 
when it is desired to predict flame velocities, only one type 
of data is available — for example, fuel-air-ratio data at 
atmospheric pressure and room temperature. It is therefore 
important to know how well the effects of initial temperature 
and oxygen concentration can be predicted from such data. 
The ethylene data are used in the following example. 

The sterie factor P calculated from the equivalence-ratio 
data and the low-temperature activation energy of 40 kilo- 
calories per gram-mole are used in the Semenov equation to 
predict flame velocities over the ranges of initial temperature 
and oxygen concentration covered by the experimental data. 
In both, cases the predicted flame velocities deviate from the 
normalized, measured flame velocities by an average of + 14 
percent, the maximum deviation being +24 percent for the 
oxygen-concentration data. (Because different methods of 
measuring flame velocity were used, differences in values for 
mixtures of the same composition and temperature were 
found. For these calculations all flame velocities are 
normalized by a simple ratio factor to the values for the 
initial-temperature data.) 

TANFOKD-PEASE EQUATION 

Tanford and Pease (refs. 11, 12, and 22) equated the 
amount of product formed in the combustion zone by a 
second-order reaction between fuel molecules and hydrogen 
atoms (or other active particles such as hydroxyl radicals or 
oxygen atoms) to the amount of product formed at the flame 
front by conversion of the fresh gas expressed in terms of 
initial conditions and flame velocity. An approximate solu- 
tion for flame velocity from this equation was given: 

(«> 

Three methods of evaluating the predictions of the Tan- 
ford and Pease equation are used herein. For all three, the 
following calculations are the same: (a) Flame temperature 
T r and active particle concentrations p t are calculated as- 
suming adiabatic thermal equilibrium by a matrix method 
(ref. 23); (b) the mean combustion-zone temperature is as- 
sumed to be 0.7 T f (ref. 11); (c) diffusion coefficients D i>m 
are calculated from D (0.7 'Tf/T 0 y- 67 , where D if 0 is 
the diffusion coefficient at initial temperature calculated by 
the Stefan-Maxwell equation (ref. 24); (d) the recombina- 
tion factor Bi is calculated by the method of Tanford (ref. 
12) for the hydrogen atoms and is assumed to be 1 for OH 
and O; (e) the ratio L m Q'n/Q. is calculated from a knowledge 
of the over-all oxidation process and the initial concentration 
of reactants. 

The three methods of evaluation differ in the calculation 
of k { . For the first method onl} r one chain carrier, the hydro- 
gen atom, is considered. For this case k H values are calcu- 


lated from single-point flame-velocity determination by 
equation (6). For the second method, H, O, and OH are 
considered to be the chain carriers and 


kiPtDi, 

Bi 


t 0 ( : 


PsDb 


B s 


-posDos, o~hpoDo, 




1.07 


( 7 ) 


Single-point k a values are calculated. For these two methods, 
the ratios of the predicted flame velocity to the experimental 
flame velocity is calculated as ( k H lk H Y 12 and ( k a /lc„y i 2 , re- 
spectively, and the average deviations of the ratios from 1 
are used as a measure of the accuracies of the predictions. 

For the third method, an Arrhenius type temperature 
dependence of the rate constant is used: 


k a =PaZ a <*P [~EJR{ 0.7 Tf)] (S) 


The following method is used to calculate an activation 
energy E a for best fit from the initial-temperature data. 
Several values of _E a are substituted in the equation and for 
each a value of ( P/PY' 2 is calculated. Then 5, the percentage 
mean deviation in ( P/P) 1 ' 2 for the group, is plotted against 
the E a used. The “fitted” value of E a is the one which 
gives the minimum 5. All other calculations by the third 
method are made by use of equations (7) and (8), measured 
flame velocities, and E a from initial-temperature data^ Tho 
accuracy is calculated as the mean deviation of {PIPY 12 
from 1. 

The predictions of the Tanford-Pease equation are evalu- 
ated by methods 1 and 2 for the maximum flame velocities 
of a series of hydrocarbons (ref. 25). Examples of the 
predictions based on k a are shown in table I. Table II 
shows that the mean deviation for 55 hydrocarbons (exclud- 
ing ethylene) is 3.3 percent considering hydrogen atoms k u 
only, and 3.5 percent considering all three active particles 
k a . These mean deviations are small enough to indicate 
either that the “rate constants” are not temperature- 
dependent, in which case the sterie factors are of the order of 
10 -3 , or that exp [E a IR(0.7 is nearly constant for the 

hydrocarbons studied, which is probable. If the activation 
energy is of the order of 10 kilocalories per gram-mole, the 
sterie factor is of the order of 10~ 2 . 

All three methods of evaluation are used with the initial- 
temperature data. The accuracy of the flame velocity 
predictions is 9 to 12 percent (table II) for the first method 
(& H ) and 6 to 8 percent for the second method. The curves 
predicted for ethylene by use of methods 2 and 3 are shown 
in figure 1. The use of method 3, which includes the 
temperature dependence of the rate constant, improves the 
prediction so that the accuracy is about 2 percent. The 
calculated activation energies and sterie factors are: for 
methane, 17 and 0.19; and for ethylene, 18 and 0.21. Tho 
calculated activation energy for ethylene is high when com- 
pared with the experimental activation energy for the 
reaction H + C 2 H 4 -s-C 2 H 5 , which is 2.6 kilocalories (calcu- 
lated from data in ref. 26). The relative predictions are not 
appreciably changed by assuming B n —\ or by assuming 
D i<m ccT m 2 which was used in reference 7 instead of 
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If llie entire range of equivalence ratios which was covered 
experimentally by the flame-velocity data is used for the 
prediction of velocities by the Tanford-Pease method, the 
accuracy of prediction is veiy poor as illustrated in figure 
2 (a) for ethylene. If the equivalence ratios are limited to a 
small range near the ratio for the maximum flame velocity 
(for example, tp=<p max ±. 0.15), the accuracy of the predictions 
is 1 to 6 percent (table I) when all active particles are con- 
sidered. However, if only hydrogen atoms are used the 
accuracy is very poor — 17 to 65 percent. 

The variations of flame velocity with oxygen concentration 
are calculated for all active particles. The accuracy of 
prediction is 4 to 6 percent. The predicted curves are 
shown in figure 3 (b). 

The ethylene data are again used as an example to show 
how well an average “rate constant” (jc a from method 2) 
from the fuel-air-ratio data predicts the effects of temperature 
and oxygen concentration on the maximum flame velocity. 
Predicted flame velocities for various temperatures differ 
from the normalized experimental values by —9 to +13 
percent with an average deviation of 6 percent. Predicted 
flame velocities for various oxygen concentrations differ 
from experimental values by +8 to +22 percent with an 
average deviation of 12 percent. 

MANSON EQUATION 

The momentum relation between the flame (or combustion 
wave) velocity and the pressure drop across a plane, 
steady-state flame front, may be expressed in the form (ref. 
27): 



Manson (ref. 13) suggested that the small pressure drop could 
be caused by the projection of hydrogen atoms into the 
unburned gas. Because the H atoms would recombine to 
H a at the unburned gas temperature, Ap was assumed to 
be one-half the equilibrium H atom pressure reduced to 
unburned gas temperature. In the present paper, Ap is cal- 



(b) Tanford and Pease square-root law (diffusion mechanism). 


Figure 3. — Concluded. Effect ol oxygen concentration on maximum 
flame velocities. Comparison of theoretically predicted curves with 
experimental results. 


culated as both | p H ^ and | (j? h +Poh jr- 

The second method of calculating A p takes into account 
three of the lightest and most abundant species which might 
be considered active particles and assumes that all of these 
would recombine in pahs. The average proportionality 
factor for a group of data is calculated by dividing the experi- 
mental flame velocity bj 5- the right side of equation (9) for 
each point and averaging the quotients. These .average 
empirical proportionality factors are designated K n for the 
first method and K s for the second. The accuracy of the 
predicted flame velocities is considered to be the average 
deviation of K B /K a or K 2 /K 2 from 1. For the maximum 
flame velocities of the different hydrocarbons, the accuracy 
of prediction by the Manson equation for either method of 
evaluation is about 4 percent (table II). For the first 
method (K s ), the proportionality factor is 0.5; that is, the 
predicted value is twice the experimental value. 

For the initial-temperature data, the predicted flame 
velocities differ from the experimental velocities by 10 to 14 
percent. When the total range of the equivalence-ratio 
data is considered, the accuracy of prediction of flame veloci- 
ties by the Manson equation and either of the methods of 
evaluating A p is poor. If the equivalence range is limited 
as for the Tanford-Pease method and with all active particles, 
predictions accurate to 2 to 4 percent are obtained (table II). 
The accuracy of prediction of flame velocities for mixtures 
containing different concentrations of oxygen varies from 
hydrocarbon to hydrocarbon. For 2,2,4-trimethylpentane, 
the accuracy is about 1 percent; for ethylene, 5 percent; and 
for propane, 15 percent. 

The ethylene data are again considered in order to deter- 
mine the accuracy with which the K 2 from the equivalence- 
ratio data predicts the effects of temperature and oxygen 
concentration on flame velocity. The deviations for the 
initial-temperature data varied from —8 to +45 percent 
with a mean deviation of 27 percent. For the oxygen data 
the deviations varied from +1 to +19 percent with a mean 
value of +10 percent. 

EMPIRICAL CORRELATIONS 

A correlation was obtained (ref. 28) by which maximum 
flame velocities of various hydrocarbons in ah could be 
predicted with an average deviation of 2 percent. The 
predicted flame velocity is calculated from the sum of the 
contributions of various H— C bonds in the fuel molecule 
according to the following relation: 

U=N a K a +N b K b +N c K c + . . . 

where N A , N B , N c , N D , N E , N F , N a , and N B are the num- 
bers of methane, primary, second aiy, tertiary, alkene, al- 
kyne, cyclohexyl, and aromatic C— H bonds, respectively, 
per unit volume of hydrocarbon-air mixture, and K A , K B , 
K c , ■ ■ ■ are the flame-speed coefficients of these bonds. 
For the special cases of C— H bonds on carbon atoms placed 
alpha to alkyne C=C bonds, a factor M=0.96 was intro- 
duced into terms representing these alpha bonds. Correla- 
tion coefficients (table III) established from 34 hydrocarbons 
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excluding ethylene then gave an average deviation of 1.9 
percent in the ratio of predicted to measured flame velocity. 

A correlation was also observed between the equivalence 
ratio for ma ximum flame velocity and the total bond dis- 
sociation energy of the fuel. As shown in table IV, the 
total bond dissociation energy per unit volume of the 
liydrocarbon-air mixture corresponding to the maximum 
flame velocitj' is nearly constant, with an average deviation of 
only 0.9 percent from the average value. This deviation is 
less than one-third of that obtained by simply assuming that 
the maximum will occur at an average equivalence ratio of 
1.15. Bond energy calculations for 37 hydrocarbons similar 
to those presented in table IV had an average deviation of 
1.6 percent. 

For the equivalence-ratio data, linear relations were found 
between flame velocity and the logarithm of the equivalence 
ratio <p, where <p<C<p ma x- An interesting variation is the plot 
of flame velocity' against the total bond dissociation energy 
of the fuel per unit volume of mixture in figure 4 for pentane, 
ethylene, and propyne mixtures with ah. It may be seen 
that the correlation is linear for the lean mixtures and that 
the line extrapolates to the lean limit for flame propagation 
(where 17=0) as determined in the 1-inch flame tube. For 
these data, determined by the tube method for three hydro- 
carbons, good estimates of the maximum flame velocities are 
obtained by reading from the craves at a total bond dis- 

TABLE III.— EMPIRICAL COEFFICIENTS FOR CALCULAT- 
ING MAXIMUM FLAME VELOCITY FROM HYDROCARBON 

STRUCTURE (REF. 28) 


Type C-E bond 


Coefficient 

Methane 

Ka 

35.2X10-W 

Primary 

Kb 

42.5 

Secondary — 

Kc 

47.5 

Tertiary 

Kd 

45.4 

Cyclohexyl 

Ka 

50.5 

Alkene 

Kb 

80.7 

Aromatic - 

K„ 

84.3 

Alkyne 

Kf 

223.9 


TABLE IV— BOND DISSOCIATION ENERGIES FOR C< AND 
C 8 HYDROCARBONS AT CONCENTRATION FOR MAXI- 
MUM FLAME VELOCITY 


Hydrocarbon 

Equiva- 
lence 
ratio 
for max- 
imum TJy 

*Pmai 

Devia- 
tion of 

*pmaz 

from 

average 

<pnax, 

percent 

Dissocia- 
tion en- 
ergy of 
hydro- 
carbon, 
kcal 

Dissocia- 
tion en- 
ergy of 
mixture 

at tptnox, 

kcal 

liter 

Devia- 

tion 

from 

average, 

percent 

g-molc 

(*) 

Butane 

1.09 

5.5 

1250 

1.93 

2.0 

Hexane 

1. 16 

0.9 

1791 

2.00 

1.5 

2-Methylpropane 

1. 11 

3.6 

1257 

1.94 

1.5 

2, 2-Dimethylbutano-- 

1.12 

2.5 

1821 

1.97 

0.0 

Butcne-I 

1. 17 

1.7 

1133 

1.98 

0.5 

He.xene-1 

1.16 

0.9 

16S0 

2.02 

2.5 

Isobutene 

1. 14 

0.9 

1147 

1.97 

0.0 

2-Mcthylpcnteno-l 

1. 19 

3.4 

1603 

1.99 

1.0 

Butyno-I 

1. 17 

1.7 

1016 

1.95 

1.0 

Hcxyne-l 

1.21 

5.0 

1564 

1.97 

0.0 

4-Methylpentyno 

1. 18 

2.5 

1574 

1.95 

1.0 

Benzene - 

1.34 

14.2 

1305 

1.96 

0.5 


1.15 

3.29 


1.96 

0.9 








“Single-bond energies Irom ref. 29. Multiple-bond energies from ref. 30: 151.2 kcal/g-mole 
for C— C in olefins and cyclic compounds; 198.5 kcal/g-mole for C=C bond. 



Figure 4. — Empirical correlation for equivalence-ratio data. 


sociation energy of 1.75 kilocalories per liter. This constant 
bond -energy value of 1.75 kilocalories per liter for this cor- 
relation does not equal the bond energy of the hydrocarbon 
concentration corresponding to the maximum flame velocity 
which was previously calculated to be 1.96 kilocalories per 
liter. Work with other flame-velocity data shows that both 
the position of the horizontal line and the slopes of the lines 
for specific hydrocarbons depend somewhat on the motliod 
of measurement of flame velocity. 

Empirical equations predicting the effect of initial mixture 
temperature on flame velocity with mean deviations of ap- 
proximately 2 percent for the range from 200° to 615° K 
are (ref. 5): 

For methane: 

£7=10+7.40X10- 5 2y- 23 


For propane: 

Z7=10-t-3.42X10 -4 2o 2, °° 


For ethylene: 

17= 10+2.59 X10- 3 To K74 


Empirical relations were found which predict the effect oJ 
oxygen concentration (for the experimental range covered] 
and limited temperature (311°-422° K) on the flame velocity 
for the three hydrocarbons studied (ref 6): 

For 2,2,4-trimethylpentane: 


17=0.133 To 1 - 40 (a — 0.120) 

For propane: 

17=0.766 To 1 ' 16 (a— 0.133) 

For ethylene: 

77=0.998 To 1 ' 18 (a-0.133) 


The accuracy of prediction was about 3 percent for ethylene 
and 2,2,4-trimethylpentane and 6 percent for propane. 


CONCLUDING REMARKS 

The semitheoretical calculations of this paper are depend' 
ent on the assumption that thermal equilibrium is attainec 
at the end of the combustion zone. Actual flame tempera' 
tures and product concentrations would differ from equilib 
rium values if (a) equilibrium is not attained, (b) any ap' 
preciable reaction takes place at the initial temperature, oi 
(c) there is a large chain branching term affecting the radica 
distribution in and ahead of the .flame zone. These calcu 
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lations are also dependent on the inherent assumptions of the 
approximate theoretical equations and the methods used to 
calculate transport properties. Within these limitations, the 
results of the calculations may be generalized as follows: 

1. The Semenov bimolecular equation (thermal mechan- 
ism), the Tanford-Pease square-root law (active-particle 
diffusion mechanism), and the Manson modification of the 
momentum-pressure-drop equation (including diffusion of 
active particles, but not chemical kinetics) will all predict 
relative changes in flame velocity caused by changes in hydro- 
carbon structure, initial temperature, equivalence ratio 
(limited range), or oxygen concentration. (The accuracies 
of the predictions are tabulated in this report.) 

2. Low-temperature activation energies may be used in 
the Semenov equation to give good relative predictions of 
flame velocities for variations in initial temperature, equiva- 
lence ratio, and oxygen concentration. Only the bimolecular 
fuel-oxygen reaction is considered in this paper; other com- 
binations such as fuel and hydrogen atoms should bo in- 
vestigated to determine whether the kinetics indicated for 
such other reactions are not more probable than for the fuel- 
oxygen case. 

3. The average “specific rate constants” k a obtained with 
the Tanford-Pease equation for molecular fuel — active par- 
ticle reactions behave as Arrhenius rate constants only with 
the data on the effect of initial temperature, not with the 
equivalence-ratio or oxj r gen-concentration data. The rela- 
tive predictions obtained by this equation are not very sensi- 
tive either to the temperature dependence assigned to the 
diffusion coefficients or to the recombination factor computed 
for hydrogen atoms; these factors do appreciably affect 
absolute predictions of flame velocity or, conversely, de- 
terminations of steric factors from experimental flame 
velocities. Regarding active particles considered, better re- 
sults are generally obtained for hydrocarbon flames when H, 
OH, and O are considered rather than H alone; this is par- 
ticularly true for studies of the effect of equivalence ratio. 

4. For the equivalence-ratio and oxygen-concentration 
data, considerably better relative predictions are obtained 
from the Manson equation when the pressure drop across the 
flame front is considered to be due to H, OH, and O than 
when only H atoms are taken into account. Further thought 
should be given as to how this pressure drop should be com- 
puted, taking into account any net diffusional flow of any 
component between the flame front and the unburned gas. 

5. The examples presented for the ethylene data show that 
an empirically determined rate constant or proportionality 
factor from equivalence-ratio data at atmospheric pressure 
and room temperature may predict the effect of initial 
temperature or oxygen concentration within approximately 
20 percent by the Semenov or Tanford-Pease equation and 
within approximately 30 percent by the Manson equation. 

6. For engineering applications, the effects of the param- 
eters studied could be estimated just as satisfactorily, and 
more easily, by one or another of the empirical correlations 
indicated, as compared with the three semitheoretical equa- 
tions considered. However, the use of the semitheoretical 
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equations in some cases reduces the number of constants 
required. 


Lewis Flight Propulsion Laboratory 
National Advisory Committee for Aeronautics 
Cleveland, Ohio, January 2, 1958 
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